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NEUTRINO PROPERTIES FROM HIGH ENERGY 
ASTROPHYSICAL NEUTRINOS 



SANDIP PAKVASA 

It is shown how high energy neutrino beams from very distant sources can be 
utihzed to learn about some properties of neutrinos such as hfetimes, mass hier- 
archy, etc. Furthermore, even mixing elements such as Ues and the CPV phase 
in the neutrino mixing matrix can be measured in principle. Pseudo-Dirac mass 
differences as small as 10~^**ey^ can be probed as well. 



1. Introduction 

We make two basic assumptions which are reasonable. The first one is that 
distant neutrino sources (e.g. AGN's and GRB's) exist; and furthermore 
with detectable fluxes at high energies (upto and beyond PeV). The second 
one is that in the not too far future, very large volume, well instrumented 
detectors of sizes of order of KM3 and beyond will exist and be operating; 
and furthermore will have (a) reasonably good energy resolution and (b) 
good angular resolution (~ 1° for muons). 

2. Neutrinos from Astrophysical Sources 

If these two assumptions are valid, then there are a number of uses these 
detectors can be put to^ . In this talk I want to focus on those that enable us 
to determine some properties of neutrinos: namely, probe neutrino lifetimes 
to W^s/eV (an improvement of 10* over current bounds), pseudo-Dirac 
mass splittings to a level of lO^^^eV^ (an improvement of a factor of 10^ 
over current bounds) and potentially even measure quantities such as Ues 
and the phase S in the MNSP matrix^. 



3. Astrophysical neutrino flavor content 

In the absence of neutrino oscillations we expect a very small v-r compo- 
nent in neutrinos from astrophysical sources. From the most discussed 
and the most likely astrophysical high energy neutrino sources^ we expect 
nearly equal numbers of particles and anti-particles, half as many v'^s as 
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u'lj^s and virtually no u'^s. This comes about simply because the neutrinos 
are thought to originate in decays of pions (and kaons) and subsequent 
decays of muons. Most astrophysical targets are fairly tenous even com- 
pared to the Earth's atmosphere, and would allow for full muon decay in 
flight. There arc some predictions for flavor independent fluxes from cos- 
mic defects and exotic objects such as evaporating black holes. Observation 
of tau neutrinos from these would have great importance. A conservative 
estimate"* shows that the prompt i^r flux is very small and the emitted flux 
is close to the ratio 1:2:0. The flux ratio of z/g : i^^ : i^^ = 1 : 2 : is 
certainly valid for those AGN models in which the neutrinos are produced 
in beam dumps of photons or protons on matter, in which mostly pion and 
kaon decay (followed by the decay of muons) supply the bulk of the neutrino 
flux. 

Depending on the amount of prompt z/— flux due to the production 
and decay of heavy flavors, there could be a small non-zero Vr component 
present. 

4. Effect of Oscillations 

The current knowledge of neutrino masses and mixings can be summarized 
as follows^. The mixing matrix elements are given to a good approximation 
with the solar mixing angle given by about 32", the atmospheric angle by 
about 45° and Ue3 < 0.17 limited by the CHOOZ bound. The mass spec- 
trum has two possibilities; normal or inverted, and with the mass differences 
given by ~ 2.10^'^ey^ and Sm2i ^ l.lO^^eV^. Since 5w?L/4,E for 

the distances to GRB's and AGN's (even for energies upto and beyond 
PeV) is very large (> 10'') the oscillations have always averaged out and 
the conversion(or survival) probability is given by 

i 

Assuming no significant matter effects enroutc. it is easy to show that the 
mixing matrix in Eq. (1) leads to a propagation matrix P, which, for any 
value of the solar mixing angle, converts a flux ratio oi : v^j, : Vt = ^ '■ 
2 : into one of 1 : 1 : 1. Hence the flavor mix expected at arrival is 
simply an equal mixture of Ve^^ix and v^- as was observed long ago"''^. If 
this universal flavor mix is confirmed by future observations, our current 
knowledge of neutrino masses and mixings is reinforced and conventional 
wisdom about the beam dump nature of the production process is confirmed 
as well. However, it would much more exciting to find deviations from it. 
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and learn something new. How can this come about? Below is a shopping 
list of a variety of ways in which this could come to pass. 

5. Deviations from Canonical Flavor Mix 

There are quite a few ways in which the flavor mix can be changed from 
the simple universal mix. 

The first and simplest is that initial flavor mix is NOT 1:2:0. This 
can happen when there arc strong magnetic fields causing muons to lose 
energy before they decay, and there exist models for neutrino production 
in AGN's in which this does happen^. In this case the v'^s have much lower 
energies compared to i^'^s and effectively the initial flavor mix is : 1 : 
and averaged out oscillations convert this into 1/2 : 1 : 1 on arrival. 

The possibility that the mass differences between neutrino mass eigen- 
states are zero in vacuum (and become non-zero only in the presence of 
matter) has been raised recently*^ . If this is true, then the final flavor mix 
should be the same as initial, namely: 1:2:0. 

Neutrino decay is another important possible way for the flavor mix to 
deviate signiflcantly from the democratic mix® .We now know that neutrinos 
have non-zero masses and non-trivial mixings, based on the evidence for 
neutrino mixings and oscillations from the data on atmospheric, solar and 
reactor neutrinos. 

If this is true, then in general, the heavier neutrinos are expected to de- 
cay into the lighter ones via flavor changing processes^". The only questions 
are (a) whether the lifetimes are short enough to be phenomenologically in- 
teresting (or are they too long?) and (b) what are the dominant decay 
modes. 

Since we are interested in decay modes which are likely to have rates (or 
lead to lifetimes) which are phenomenologically interesting, we can rule out 
several classes of decay modes immediately. For example, the very strong 
constraints on radiative decay modes and on three body modes such as 
v render them as being uninteresting. 

The only decay modes which can have interestingly fast decays rates are 
two body modes such as i/i Vj + x where x is a very light or massless par- 
ticle, e.g. a Majoron. In general, the Majoron is a mixture of the Gelmini- 
Roncadelli^^ and Chikasige-Mohapatra-Peccei^^ type Majorons. The effec- 
tive interaction is of the form: 



i/^(a -I- h^5)va X 



(2) 
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giving rise to decay: 



■a 



Uf) {or yp) + X 



(3) 



where x is a massless, spinless particle; Ua and Ujs are mass eigenstates 
which may be mixtures of flavor and sterile neutrinos. Explicit models of 
this kind which can give rise to fast neutrino decays have been discussed^'^. 
These models are unconstrained by /z and r decays which do not arise 
due to the AL = 2 nature of the coupling. The couplings of and are 
constrained by the limits on multi-body tt, K decays, and on /i — e university 
violation in tt and K decays^**, but these bounds allow fast neutrino decays. 

There are very interesting cosmological implications of such couplings. 
The details depend on the spectrum of neutrinos and the scalars in the 
model. For example, if all the neutrinos are heavier than the scalar; the 
relic neutrino density vanishes today, and the neutrino mass bounds from 
CMB and large scale structure are no longer operative, whereas future 
measurements in the laboratory might flnd a non-zero result for a neutrino 
mass . If the scalars are heavier than the neutrinos, there are signatures 
such as shifts of the nth multipole peak (for large n) in the CMB There 
are other implications as well, such as the number of relativistic degrees of 
freedom(or effective number of neutrinos) being different at the BEN and 
the CMB eras. The additional degrees of freedom should be detectable in 
future CMB measurements. 

Direct limits on such decay modes are also very weak. Current bounds 
on such decay modes are as follows. For the mass eigenstate ux, the limit 
is about 



based on observation of Vf.s from SN1987A (assuming CPT invariance). 
For z/2, strong limits can be deduced from the non-observation of solar 
anti- neutrinos in KamLAND"'^*, a more general bound is obtained from an 
analysis of solar neutrino data^*^ leads to a bound given by: 



For 1^3, in case of normal hierarchy, one can derive a bound from the atmo- 
spheric neutrino observations of upcoming neutrinos'^" : 



n > 10^ sec/eV 



(4) 



T2 > 10~'^ sec/eV 



(5) 



T3 > 10-^° sec/eV 



(6) 



The strongest lifetime limit is thus too weak to eliminate the possibility 
of astrophysical neutrino decay by a factor about 10'' x (L/lOO Mpc) x(10 
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TeV/E). It was noted that the disappearance of all states except i/i would 
prepare a beam that conld in principle be used to measure elements of the 
neutrino mixing matrix, namely the ratios C/g]^ : U^i : U^i^^. The possibility 
of measuring neutrino lifetimes over long baselines was mentioned in Ref.^^, 
and some predictions for decay in four- neutrino models were given in Ref.^^. 
The particular values and small uncertainties on the neutrino mixing pa- 
rameters allow for the first time very distinctive signatures of the effects of 
neutrino decay on the detected flavor ratios. The expected increase in neu- 
trino lifetime sensitivity (and corresponding anomalous neutrino couplings) 
by several orders of magnitude makes for a very interesting test of physics 
beyond the Standard Model; a discovery would mean physics much more 
exotic than neutrino mass and mixing alone. Neutrino decay because of 
its unique signature cannot be mimicked by either different neutrino flavor 
ratios at the source or other non-standard neutrino interactions. 

A characteristic feature of decay is its strong energy dependence: 
exp(— Lm/i?T), where t is the rest-frame lifetime. For simplicity, con- 
sider the case that decays are always complete, i.e., that these exponential 
factors vanish. 

The simplest case (and the most generic expectation) is a normal hi- 
erarchy in which both and 1^2 decay, leaving only the lightest stable 
eigenstate vi. In this case the flavor ratio is U^^ : U'^i : U^i^^- Thus if 

Ue3 = 

(t>ve : t/);.^ : 01/, ^ 5 : 1 : 1, (7) 

where we used the neutrino mixing parameters given above^. Note that 
this is an extreme deviation of the flavor ratio from that in the absence of 
decays. It is difficult to imagine other mechanisms that would lead to such 
a high ratio of i^p to v^^ . In the case of inverted hierarchy, 1/3 is the lightest 
and hence stable state, and so^ 

: K ■■ 0-. = : (7^3 : f/^g = : 1 : 1. (8) 

If Uez = and Oatm = 45°, each mass eigenstate has equal and Vr com- 
ponents. Therefore, decay cannot break the equality between the (f)^,^ and 
(/>y, fluxes and thus the 0i/^ : 0^,^ ratio contains all the useful information. 
The effect of a non-zero Ues on the no-decay case of 1 : 1 : 1 is negligible. 

When Ue3 is not zero, and the hierarchy is normal, it is possible to obtain 
information on the values of Ue3 as well as the CPV phase (5^^. The flavor 
ratio e/ii varies from 5 to 15 (as Ues goes from to 0.2) for cos6 = +1 but 
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from 5 to 3 for cosS = — 1. The ratio t//U varies from 1 to 5 (cos 5 = +1) 
or 1 to 0.2 (cos (5 — —1) for the same range of 1/^3- 

If the decays are not complete and if the daughter does not carry the 
full energy of the parent neutrino; the resulting flavor mix is somewhat 
different but any case it is still quite distinct from the simple 1:1:1 mix^. 

If the path of neutrinos takes them thru regions with significant mag- 
netic fields and the neutrino magnetic moments are large enough, the flavor 
mix can be affected^^. The main effect of the passage thru magnetic field 
is the conversion of a given hclicity into an equal mixture of both helicity 
states. This is also true in passage thru random magnetic fields^^. 

If the neutrino are Dirac particles, and all magnetic moments are com- 
parable, then the effect of the spin-flip is to simply reduce the overall flux 
of all flavors by half, the other half becoming the sterile Dirac partners. If 
the neutrinos are Majorana particles, the flavor composition remains 1 : 1 
: 1 when it starts from 1 : 1 : 1, and the absolute flux remains unchanged. 

What happens when large magnetic fields are present in or near the neu- 
trino production region? In case of Dirac neutrinos, there is no difference 
and the outcoming flavor ratio remains 1 : 1 : 1, with the absolute fluxes 
reduced by half. In case of Majorana neutrinos, since the initial flavor mix 
is no longer universal but is z^e : : ^^-r ~ 1 : 2 : 0, this is modified but it 
turns out that the final(post-oscillation) flavor mix is still 1:1:1! 

Other neutrino properties can also affect the neutrino flavor mix and 
modify it from the canonical 1 : 1 : 1. If neutrinos have flavor(and equiva- 
lence principle) violating couplings to gravity (FVG), or Lorentz invariance 
violating couplings; then there can be resonance effects which make for one 
way transitions(analogues of MSW transitions) e.g. — > Vt but not vice 
versa^^'^^. In case of FVG for example, this can give rise to an anisotropic 
deviation of the Vfi/vT ratio from 1, becoming less than 1 for events coming 
from the direction towards the Great Attractor, while remaining 1 in other 
directions'^. 

Another possibility that can give rise to deviations of the flavor mix from 
the canonical 1 : 1 : 1 is the idea of neutrinos of varying mass(MaVaNs). In 
this proposal'^, by having the dark energy and neutrinos(a sterile one to be 
specific) couple, and track each other; it is possible to relate the small scale 
2 X 10~^ eV required for the dark energy to the small neutrino mass, and 
furthermore the neutrino mass depends inversely on neutrino density, and 
hence on the epoch. As a result, if this sterile neutrino mixes with a flavor 
neutrino, the mass difference varies along the path, with potential resonance 
enhancement of the transition probability into the sterile neutrino, and 
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thus change the flavor mix^°. For example, if only one resonance is crossed 
enroutc, it can lead to a conversion of the lightest (mostly) flavor state 
into the (mostly) sterile state, thus changing the flavor mix to 1 — U^i : 
1 — U^-^ : 1 — U^-^ « 1/3 : 1 : 1, in case of inverted hierarchy and similarly 
pa 2 : 1 : 1 in case of normal hierarchy. 

Complete quantum decoherence would give rise to a flavor mix given 
by 1 : 1 : 1, which is identical to the case of averaged out oscillations as 
we saw above. The distinction is that complete decoherence always leads 
to this result; whereas averaged out oscillations lead to this result only in 
the special case of the initial flavor mix being 1 : 2 : 0. To find evidence for 
decoherence, therefore, requires a source which has a different flavor mix 
. One possible practical example is a source which emits i^'^s by decay of 
neutrons, and hence no f'^s at all, with an initial flavor mix of 1 : : 0. 
In this case decoherence leads to the universal 1:1:1 mix whereas the 
averaged out oscillations lead to 3 : 1 : 1^^. The two cases can be easily 
distinguished from each other. 

If each of the three neutrino mass eigenstates is actually a doublet with 
very small mass difference (smaller than 10~^eV), then there are no current 
experiments that could have detected this. Such a possibility was raised 
long ago"*^. It turns out that the only way to detect such small mass 
differences {W^^^eV^ > 5m^ > 10~^^eV^) is by measuring flavor mixes of 
the high energy neutrinos from cosmic sources. Relic supernova neutrino 
signals and AGN neutrinos are sensitive to mass difference squared down 

to 10-2061^2 33 

Let {v^ ,1^2 , ^t'^^i T^z) denote the six mass eigenstates where u'^ 
and v~ are a nearly degenerate pair. A 6x6 mixing matrix rotates the mass 
basis into the flavor basis. In general, for six Majorana neutrinos, there 
would be fifteen rotation angles and fifteen phases. However, for pseudo- 
Dirac neutrinos, Kobayashi and Lim^^ have given an elegant proof that the 
6x6 matrix Vkl takes the very simple form (to lowest order in Sw? /w?: 

''-=(o^J-(:^:-r4)' 

where the 3x3 matrix U is just the usual mixing matrix determined by the 
atmospheric and solar observations, the 3x3 matrix Ur is an unknown uni- 
tary matrix and Vi and V2 are the diagonal matrices Vi = diag (1, 1, l)/\/2, 

and V2=diag(e^*'*% e^*"^^, e~*'^^)/A/2, with the (fn being arbitrary phases. 
As a result, the three active neutrino states are described in terms of 
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VaL = Uccj + ii^j ) . 

The flavors deviate from the democratic value of i by 



5Pe 
5Pr 



4X1 + 4X2 



1 3 1 
gXl + g X2 + 2 X3 



(10) 



where Xi = sin^((5mfZ//4£').The flavor ratios deviate from 1:1:1 when one 
or two of the pseudo-Dirac oscillation modes is accessible. In the ultimate 
limit where L/E is so large that all three oscillating factors have averaged 
to ^, the flavor ratios return to 1 : 1 : 1, with only a net suppression of the 
measurable flux, by a factor of 1/2. 



6. Cosmology w^ith Neutrinos 

If the oscillation phases can indeed be measured for the very small mass 
differences by the deviations of the flavor mix from 1 : 1 : 1 as discussed 
above, the following possibility is raised. It is a fascinating fact that non- 
averaged oscillation phases, S(f)j = Srrijt/Ap, and hence the factors Xj, are 
rich in cosmological information^^ Integrating the phase backwards in 
propagation time, with the momentum blue-shifted, one obtains 



= / 

Jo 



4po 



(12) 



where I is given by 



-L 



l+Ze 



Ze is the red-shift of the emitting source, and H^^ is the Hubble time, 
known to 10% . This result holds for a flat universe, where Q.m + Ha = 1, 
with and JIa the matter and vacuum energy densities in units of the 
critical density. The integral I is the fraction of the Hubble time available 
for neutrino transit. For the presently preferred values = 0.3 and 
Qa = 0.7, the asymptotic {ze oo) value of the integral is 0.53. This limit 
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is approached rapidly: at Ze = 1 (2) the integral is 77% (91%) saturated. 
For cosmologically distant {z^ > 1) sources such as gamma-ray bursts, non- 
averaged oscillation data would, in principle, allow one to deduce 6m^ to 
about 20%, without even knowing the source red-shifts. Known values of 
Qm and J^a might allow one to infer the source rcdshifts Zg, or vice- versa. 

This would be the first measurement of a cosmological parameter with 
particles other than photons. An advantage of measuring cosmological 
parameters with neutrinos is the fact that flavor mixing is a microscopic 
phenomena and hence presumably free of ambiguities such as source evo- 
lution or standard candle assumptions^^'"^''. Another method of measuring 
cosmological parameters with neutrinos is given in Ref.^^. 

7. Experimental Flavor Identification 

It is obvious from the above discussion that flavor identification is crucial for 
the purpose at hand. In a water cerenkov detector flavors can be identified 
as follows. 

The fiux can be measured by the fi's produced by the charged current 
interactions and the resulting fi tracks in the detector which are long at 

these energies, i/'s produce showers by both CC and NC interactions. The 
total rate for showers includes those produced by NC interactions of i/^s 
and v'^s as well and those have to be (and can be) subtracted off to get the 
real fiux of v'^s. Double-bang and lollipop events are signatures unique to 
tau neutrinos, made possible by the fact that tau leptons decay before they 
lose a significant fraction of their energy. A double-bang event consists of a 
hadronic shower initiated by a charged-current interaction of the Vr followed 
by a second energetic shower from the decay of the resulting tau lepton*. 
A lollipop event consists of the second of the double-bang showers along 
with the reconstructed tau lepton track (the first bang may be detected or 
not). In principle, with a sufficient number of events, a fairly good estimate 
of the flavor ratio : v^^ : can be reconstructed, as has been discussed 
recently. Deviations of the flavor ratios from 1:1:1 due to possible 
decays are so extreme that they should be readily identifiable^®. Future 
high energy neutrino telescopes, such as Icccubc^'^, will not have perfect 
ability to separately measure the neutrino flux in each flavor. However, the 
situation is salvagable. In the limit of — symmetry the fluxes for 
and Vt are always in the ratio 1:1, with or without decay. This is useful 
since the v-r flux is the hardest to measure. 

Even when the tau events are not at all identifiable, the relative number 
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of shower events to track events can be related to the most interesting 
quantity for testing decay scenarios, i.e., the v,, to ratio. The precision of 
the upcoming experiments should be good enough to test the extreme flavor 
ratios produced by decays. If electromagnetic and hadronic showers can be 
separated, then the precision will be even better"^^. Comparing, for example, 
the standard flavor ratios of 1 : 1 : 1 to the possible 5:1:1 generated by 
decay, the more numerous electron neutrino flux will result in a substantial 
increase in the relative number of shower events. The measurement will be 
limited by the energy resolution of the detector and the ability to reduce 
the atmospheric neutrino background which drops rapidly with energy and 
should be negligibly small at and above the PeV scale. 

8. Discussion and Conclusions 

The flux ratios we discuss are energy-independent because we have assumed 
that the ratios at production are energy-independent, that all oscillations 
are averaged out, and that all possible decays are complete. In the standard 
scenario with only oscillations, the flnal flux ratios are (i>y^_ : 4>v,, '■ 'Pv^ = 
1:1:1. In the cases with decay, we have found rather different possible 
flux ratios, for example 5 : 1 : 1 in the normal hierarchy and : 1 : 1 in 
the inverted hierarchy. These deviations from 1:1:1 are so extreme that 
they should be readily measurable. 

If we are very fortunate^^, we may be able to observe a reasonable 
number of events from several sources (of known distance) and/or over a 
sufflcient range in energy. Then the resulting dependence of the flux ratio 
[ve/Vfj) on L/E as it evolves from say 5 (or 0) to 1, can be clear evidence of 
decay and further can pin down the actual lifetime instead of just placing 
a bound. 

To summarize, we suggest that if future measurements of the flavor mix 
at earth of high energy astrophysical neutrinos find it to be 

<l>.Jct>,Jct>,^=a/l/l; (13) 

then 

(i) awl (the most boring case) confirms our knowledge of the MNSP^ 
matrix and our prejudice about the production mechanism; 

(ii) a Ri 1/2 indicates that the source emits pure i/^s and the mixing 
is conventional; 

(iii) a « 3 from a unique direction, e.g. the Cygnus region, would be 
evidence in favour of a pure production as has been suggested 
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recently 

(iv) a > 1 indicates that neutrinos are decaying with normal hierarchy; 
and 

(v) a would mean that neutrino decays are occuring with inverted 

hierarchy; 

(vi) Values of a which cover a broader range (3 to 15) and deviation 
of the /x/r ratio from 1 (between 0.2 to 5) can yield valuable in- 
formation about Ue3 and cos S. Deviations of a which are less ex- 
treme(between 0.7 and 1.5) can also probe very small pseudo-Dirac 

(5m^ (smaller than lO^^^eV^). 

Incidentally, in the last three cases, the results have absolutely no depen- 
dence on the initial flavor mix, and so are completely free of any dependence 
on the production model. So either one learns about the production mecha- 
nism and the initial flavor mix, as in the first three cases, or one learns only 
about the neutrino properties, as in the last three cases. To summarise, 
the measurement of neutrino flavor mix at neutrino telescopes is absolutely 
essential to uncover new and interesting physics of neutrinos. In any case, 
it should be evident that the construction of very large neutrino detectors 
is a "no lose" proposition. 
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